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Abstract 

This  research  proposes  a  model  that  predicts  the  effect  of  the  anode  diffusion  layer  and  membrane  properties  on  the  electrochemical  performance 
and  methanol  crossover  of  a  direct  methanol  fuel  cell  (DMFC)  membrane  electrode  assembly  (MEA).  It  is  an  easily  extensible,  lumped  DMFC 
model.  Parameters  used  in  this  design  model  are  experimentally  obtainable,  and  some  of  the  parameters  are  indicative  of  material  characteristics. 
The  quantification  of  these  material  parameters  builds  up  a  material  database.  Model  parameters  for  various  membranes  and  diffusion  layers  are 
determined  by  using  various  techniques  such  as  polarization,  mass  balance,  electrochemical  impedance  spectroscopy  (EIS),  and  interpretation  of 
the  response  of  the  cell  to  step  changes  in  current.  Since  the  investigation  techniques  cover  different  response  times  of  the  DMFC,  processes  in 
the  cell  such  as  transport,  reaction  and  charge  processes  can  be  investigated  separately.  Properties  of  single  layers  of  the  MEA  are  systematically 
varied,  and  subsequent  analysis  enables  identification  of  the  influence  of  the  layer’s  properties  on  the  electrochemical  performance  and  methanol 
crossover.  Finally,  a  case  study  indicates  that  the  use  of  a  membrane  with  lower  methanol  diffusivity  and  a  thicker  anode  micro-porous  layer  (MPL) 
yields  MEAs  with  lower  methanol  crossover  but  similar  power  density. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  is  seen  as  a  promis¬ 
ing  power  source  for  mobile  electronic  applications.  However, 
major  problems  that  cause  low  performance  of  the  DMFC  like 
high  methanol  crossover  and  low-methanol  oxidation  activity 
have  not  yet  been  solved.  To  improve  DMFC  performance,  many 
experimental  studies  on  the  variation  of  material  and  layers 
inside  the  DMFC  have  been  carried  out  [1-7].  These  investi¬ 
gations  have  been  matched  by  several  works  on  modelling  that 
mostly  focused  on  the  effect  of  specific  materials  and  operating 
conditions  on  the  DMFC.  A  recent  review  of  DMFC  modelling 
[8]  categorises  the  model-based  research. 

A  frequent  approach  to  research  on  membrane  electrode 
assembly  (MEA)  design  concentrates  on  varying  one  layer  and 
neglecting  the  effects  of  other  DMFC  layers  on  its  perfor¬ 
mance.  For  a  systematic  approach  to  finding  an  optimum  set 
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of  MEA  materials,  the  screening  of  a  large  number  of  combina¬ 
tions  of  materials  or  layers  would  be  required.  Conducting  such 
investigations  experimentally  is  expensive  and  time-consuming. 
Therefore,  a  model-based  performance  prediction  and  optimiza¬ 
tion  are  of  great  interest  for  the  MEA  designer. 

This  research  proposes  an  easily  extendable  DMFC  MEA 
design  model.  Three  main  points  lead  to  the  MEA  design 
model.  First,  a  model  is  designed  that  is  able  to  reproduce  the 
experimental  data.  Parameters  used  in  this  design  model  are 
experimentally  obtainable,  and  some  of  these  parameters  are 
indicative  of  material  characteristics.  The  quantification  of  these 
material  parameters  builds  up  a  material  database.  Second,  pro¬ 
cesses  in  the  DMFC  such  as  transport,  reaction  and  charge  are 
separated,  allowing  their  influence  on  the  DMFC  behaviour  to 
be  independently  studied.  This  is  undertaken  by  applying  dif¬ 
ferent  investigative  methods  such  as  polarization,  current-steps 
and  electrochemical  impedance  spectroscopy  (EIS),  rather  than 
using  just  a  steady-state  method  (polarization)  which  contains 
the  sum  of  all  effects.  Finally,  the  MEA  layer  properties  are  var¬ 
ied  and  key  parameters  for  reproducing  the  features  of  each  layer 
are  identified. 
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Nomenclature 

As  geometric  electrode  area  of  DMFC  (m2) 

AWad  areal  weight  of  anode  diffusion  layer  (g  m-2) 
cpt  active  surface  concentration  of  Pt 
(=0.117  molm-2) 

cru  active  surface  concentration  of  Ru 
(=0.165  molm-2) 

cch3oh  methanol  concentration  inside  anode  compart¬ 
ment  (mol  m-3) 

cch3oh  methanol  concentration  of  inlet  feed  (mol  m-3) 
cch3oh  methanol  concentration  in  anode  catalyst  layer 
(mol  m-3) 

cc%OH  methanol  concentration  at  interface  of  anode 
diffusion  layer  and  anode  micro-porous  layer 
(molm-3) 

CAC  anode  double  layer  capacitance  (=3348  F  m-2) 
Ccc  cathode  double  layer  capacitance  (=907  F  m-2) 
dAD  diffusion  layer  thickness  excluding  MPL  (m) 
dA^  effective  diffusion  layer  thickness  excluding  MPL 

(m) 

dAMPL  thickness  of  anode  micro-porous  layer  (m) 

^hyd  thickness  of  GDL  governed  by  convective  flow 

(m) 

dM  thickness  of  membrane  (m) 

A:h3oh  diffusion  coefficient  of  methanol  in  water  at 
333  K  (=3.79  x  10-9m2s-1) 

^ch3  oh  diffusion  coefficient  of  methanol  in  membrane  at 
333  K  (m2  s-1) 

F  Faraday  constant  (=96485  C  mol- 1 ) 

Fa  flow  rate  of  flow  entering  anode  compartment 
(m3  s-1) 

geo;  (&oh)  inhomogeneity/interaction  factor  for 
Frumkin/Temkin  adsorption  on  Pt  =  1 1  and 
Ru  =  0.43 

/cell  total  cell  current  density  (A  m-2) 
kAD  effective  mass  transport  coefficient  in  anode  dif¬ 
fusion  layer  (ms-1) 

^ampl  effective  mass  transport  coefficient  in  anode 
micro-porous  layer  (ms-1) 

r\i  reaction  rate  for  anode  reaction  step  i 
(mol  m  2  s-1) 

taio  reaction  rate  constant  for  anode  reaction  step  1 
(=1.6  x  10-4ms-1) 

fA20  reaction  rate  constant  for  anode  reaction  step  2 
(=7.2  x  10-4molm-2  s-1) 

Luo  reaction  rate  constant  for  backward  reaction  of 
anode  reaction  step  2  (=9.91  x  104  molm-2  s-1) 
ta30  reaction  rate  constant  for  anode  reaction  step  3 
(=0. 19  mol  m  2  s-1) 

rc  cathode  reaction  rate  (mol  m-2  s  - 1 ) 
rco  cathode  reaction  rate  constant 
(=4  x  10-6molm-2  s-1) 

R  universal  gas  constant  (=8.314  J  mol-1  K-1) 

Re i  ohmic  resistance  of  membrane  (Q  m2) 


t 

time  (s) 

T 

temperature  of  DMFC  (=343  K) 

TjO 

uo 

standard  cell  voltage  (=1.213  V) 

fTcell 

cell  voltage  (V) 

VA 

volume  of  anode  compartment  (m3) 

yAC 

volume  of  anode  catalyst  layer  (m3) 

AMPL 

^carbon 

carbon  loading  of  micro-porous  layer  (kgm-2) 

^carbon  eff  effective  carbon  loading  of  micro-porous  layer 

VAD 

XpxFE 

(kgm  2) 

PTFE  content  in  anode  diffusion  layer  (excluding 
MPL) 

AMPL 

XpTFE 

PTFE  content  in  micro-porous  layer 

Greek  symbols 

«A2  (o?c)  charge  transfer  coefficient  for  anode  reaction 

step  2  (cathode  reaction),  =0.5  (=0.18) 

Pco  (Pou)  symmetry  parameter  for  Frumkin/Temkin 

gAO 

adsorption  on  Pt  (Ru)  (=0.5) 

porosity  of  anode  diffusion  layer  excluding 

micro-porous  layer 

gAMPL 

porosity  of  micro-porous  layer 

'?A 

anode  overpotential  (V) 

nc 

cathode  overpotential  (V) 

&CO 

surface  coverage  of  Pt  with  COads 

#OH 

surface  coverage  of  Ru  with  OHacis 

Pcarbon 

bulk  density  of  carbon  (=1.7  x  103  kgm-3) 

Pptfe 

bulk  density  of  PTFE  (=2. 19  x  103  kg  m-3) 

rAMPL 

tortuosity  of  micro-porous  layer 

2.  Experimental 

2.7.  ME  A  preparation 

In-house  MEAs  used  either  Nation  115  (DuPont)  or  hydro¬ 
carbon  membranes.  The  catalyst  was  Pt-Ru  (1:1)  black 
(HiSpec  6000,  Johnson  Matthey)  and  Pt  black  (HiSpec  1000, 
Johnson  Matthey)  for  anode  and  cathode,  respectively.  Vari¬ 
ous  anode  gas-diffusion  layers  (GDLs)  (SGL  carbon  group) 
were  used:  31BC,  31DC,  31BA  and  35BC.  They  differed 
in  poly(tetrafluoroethylene)  (PTFE)  content  and  micro-porous 
layer  (MPL).  For  the  cathode,  10DA  or  35DC  (SGL  carbon 
group)  were  used;  both  GDLs  contain  20wt.%  PTFE.  On  the 
surface  of  the  10DA  cathode  diffusion  layer,  a  1 : 1  (w/w)  mixture 
of  carbon  (Vulcan  XC-72,  Cabot)  and  PTFE  (DuPont)  dispersed 
in  isopropyl  alcohol  was  applied  in  order  to  form  a  micro- 
porous  carbon  layer  (MPL).  The  carbon  loading  of  the  MPL  was 
controlled  to  be  1.3  mg  cm-2.  Catalyst  inks,  which  consisted 
of  appropriate  amounts  of  unsupported  catalysts,  Nation  115 
solution  and  isopropyl  alcohol,  were  homogenized.  The  anode 
catalyst  layer  was  prepared  by  spraying  the  ink  on  to  the  mem¬ 
brane.  The  Pt-Ru  loading  was  8  mg  cm-2,  and  the  weight  ratio 
of  Nation  115  to  catalyst  was  0.15.  The  cathode  catalyst  layer 
was  prepared  by  spraying  Pt-black  ink  on  to  the  membrane.  The 
layer  had  a  mean  Pt  loading  of  6  mg  cm-2  and  the  Nation  115 
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to  catalyst  weight  ratio  was  0.16.  The  MEA  was  assembled  by 
hot  pressing  the  parts  at  398  K  and  500  kgf  cm-2. 

The  cell  active  area  was  either  10  x  10  4m2  or 
26xl0-4m2.  The  MEA  was  sandwiched  between  two 
plates  with  serpentine  flow  channels.  Before  recording  any 
data,  the  cell  was  activated  for  several  days. 

2.2.  Electrochemical  characterization 

All  measurements  were  conducted  with  a  Fuel  Cell  Test  Sta¬ 
tion  (Wonatech,  Korea).  For  impedance  measurements,  an  IM6 
from  Zahner  Elektrik  was  used.  Electrical  heaters  and  a  thermo¬ 
couple  were  embedded  in  the  bipolar  plates  for  controlling  the 
desired  operating  temperature  of  343  K.  A  pump  fed  the  anode 
with  methanol  solution. 

On  the  anode,  stoichiometric  feed  of  methanol  was  provided. 
However,  for  current  densities  below  1000  Am-2  a  minimum 
flow  rate  was  applied;  the  minimum  flow  rate  was  adjusted  to 
be  identical  to  the  stoichiometric  flow  rate  at  1000  Am-2. 

During  polarization,  current  step  and  mass  balance  experi¬ 
ments,  the  cathode  was  fed  with  air  at  atmospheric  pressure  at 
a  constant  flow  rate  of  500  x  10-6  m3  min-1  (298  K,  10-5  Pa). 
During  impedance  measurements,  hydrogen  at  a  flow  rate  of 
100  x  10-6  m3  min-1  (298  K,  10-5  Pa)  was  fed  to  the  cathode. 
During  recording  of  the  polarization  curve,  the  current  den¬ 
sity  was  increased  by  100 Am-2  every  70s  for  low  currents 
(<2  A);  for  higher  currents,  the  current  density  was  increased 
every  30s  by  200 Am-2.  Except  for  measurements  close  to 
the  open-circuit  voltage  (OCV)  at  which  the  interaction  of  the 
adsorbed  methanol  intermediates  may  prevent  a  rapid  establish¬ 
ment  of  a  steady-state  [9],  the  cell  achieved  a  steady-state  voltage 
within  10  s.  For  investigation  of  the  cell  voltage  response  to  step 
changes  in  current,  three  different  step  changes  were  applied. 
Two  were  of  equal  step  size  with  the  first  one  in  the  low  current 
density  regime  (1800  —>  1000  A  m-2)  and  the  second  in  the  high 
current  density  regime  (3800  — >  3000  Am-2).  The  third  exper¬ 
iment  was  conducted  in  the  low  current  density  regime  using  a 
smaller  step  size  (1800  —>  1000  A  m-2).  Each  current  was  estab¬ 
lished  for  3  min  to  guarantee  that  no  interference  between  the 
single  measurements  would  take  place.  The  current  step  exper¬ 
iments  were  repeated  three  times.  EIS  were  recorded  at  1000, 
2000  and  3000  Am-2  in  a  frequency  range  from  100  to  0.1  Hz. 
Measurements  were  performed  under  galvanostatic  control  of 
the  cell.  The  amplitude  of  the  sinusoidal  voltage  signal  was 
5  mV.  The  cathode  acted  as  a  reversible  hydrogen  electrode  and 
the  anode  spectra  were  measured  as  a  two-electrode  set-up  using 
the  anode  and  the  reversible  hydrogen  electrode.  Mass  balances 
were  conducted  to  determine  the  methanol  loss  due  to  crossover. 
Balancing  was  done  by  measuring  the  mass  change  and  concen¬ 
tration  of  the  anode  inlet  and  outlet  flows,  as  well  as  the  mass 
of  water  condensed  at  the  cathode  outlet,  at  which  the  flow  was 
cooled  to  273  K.  The  methanol  concentration  was  determined 
by  refractive  index  measurements  using  an  Atago  refractome- 
ter  (RX5000a).  The  measurements  were  conducted  over  a  time 
interval  of  1 80  min  at  an  anode  and  cathode  flow  rate  of  3  x  10-6 
and  500  x  10-6  m3  min-1  (298  K,  10-5  Pa),  respectively.  Since 
methanol  is  volatile,  this  method  may  give  results  of  medium 
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Fig.  1 .  Graphical  representation  of  phenomena  covered  by  DMFC  model  (terms 
in  brackets  are  not  taken  into  consideration). 

accuracy.  Alternative  methods  using  CO2  measurements  for 
determination  of  crossover  are  given  in  [10,11]. 

3.  Model 

3.1.  Model  set-up 

The  DMFC  model  is  a  non-linear  mathematical  model  that 
accounts  for  the  following  phenomena  (see  Fig.  1): 

•  diffusive  mass  transport  of  methanol  through  the  anode  dif¬ 
fusion  layer, 

•  oxidation  of  methanol  in  the  anode  catalyst  layer, 

•  formation  of  the  intermediates  CO  (adsorbed  on  Pt)  and  OH 
(adsorbed  on  Ru)  in  the  anode  catalyst  layer, 

•  electrochemical  reduction  of  oxygen  at  the  cathode  catalyst 
layer, 

•  methanol  crossover, 

•  undesired  electrochemical  oxidation  of  methanol  at  the  cath¬ 
ode  catalyst  layer. 

The  model  includes  the  following  assumptions: 

•  ohmic  drops  in  the  current-collectors  and  electric  connections 
are  negligible, 

•  the  fuel  cell  is  operated  isothermally, 

•  oxygen  and  carbon  dioxide  do  not  diffuse  into  the  polymer 
electrolyte  membrane  (PEM), 

•  oxygen  is  fed  in  excess,  i.e.,  oxygen  conversion  in  the  cathode 
compartment  is  negligible, 

•  the  anode  side  water  concentration  is  constant  (excess  com¬ 
ponent  in  a  liquid  mixture), 

•  on  the  anode  side,  a  pure  liquid  phase  mixture  is  assumed; 
CO2  is  postulated  to  dissolve  instantaneously, 

•  the  anode  compartment  is  perfectly  mixed, 

•  mass  transport  within  the  diffusion  layers  is  fast,  i.e.,  in  a 
quasi-steady  state, 

•  mass  transport  resistance  in  the  catalyst  layer  is  negligible, 

•  cathode  side  methanol  reacts  instantaneously  with  oxygen, 
i.e.,  its  concentration  on  the  cathode  is  zero, 
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•  no  pressure-driven  transport  through  the  membrane, 

•  electro-osmotic  drag  of  methanol  is  negligible  compared  with 
methanol  diffusive  transport, 

•  mass  transport  in  the  membrane  is  in  a  quasi-steady  state  [9], 

•  a  three- step  reaction  mechanism  is  used  to  describe  the  kinet¬ 
ics  of  the  electrochemical  oxidation  of  methanol  [12]. 


The  following  set  of  equations  holds  for  MEAs  without  an 
anode  MPL: 


dCCH3OH 
d  t 


Z7  A 

^  *  A, in 
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RT 
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RT 
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r A3  =  C A30  exp  [(1  -  fico)gco(0co  —  0.5)]^co^oh  (9) 


rc  =  -^co  exp 


(1  -  ac)F 
RT 


0c 


(10) 


UCQ n  =  Uq  -  r)A  +  0C  -  ^eHcell 


(11) 


calculation  of  the  overall  cell  voltage  (/ce li-  A  detailed  variable 
and  parameter  notation  is  given  in  the  nomenclature.  In  addition, 
Fig.  1  illustrates  the  phenomena  covered  by  this  set  of  equations 
and  their  interaction.  If  the  diffusion  layer  contains  a  MPL,  an 
additional  description  for  this  layer  is  required,  since  its  prop¬ 
erties  differ  strongly  from  that  of  the  main  GDL.  As  with  the 
GDL,  a  quasi-steady  state  concentration  profile  is  assumed  in 
the  MPL.  The  concentration  at  the  interface  between  the  GDL 
and  MPL  is  denoted  as  c^qH,  and  the  mass  transfer  coeffi¬ 
cient  for  the  MPL  as  kAMPL.  Eqs.  (1)  and  (2)  are  replaced  by  the 
following: 


dCCH3OH 
d  t 


=  77a  (T 
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The  model  should  reproduce  the  polarization,  EIS  and 
current- step  response  experiments  with  one  set  of  equations 
and  parameters.  All  three-investigation  methods  use  a  change 
in  cell  current  as  an  input  for  the  experiment  (Fig.  2).  For 
steady- state  investigations,  the  current  is  changed  step-wise  and 
only  the  respective  steady- state  voltage  is  recorded;  the  cur¬ 
rent  step  investigation  evaluates  the  transient  voltage  response 
to  a  current  step.  Finally,  in  impedance  spectroscopy,  the 
current  is  varied  sinusoidally,  and  the  ratio  between  the  volt¬ 
age  and  current,  i.e.,  the  impedance,  is  recorded.  For  the 
steady- state  and  current  step  investigations,  the  set  of  equa¬ 
tions  can  be  used  in  its  original  form  (Eqs.  (1)— (14)).  To 
model  the  recorded  anode  impedance  spectra,  the  anode  bal¬ 
ance  equations  (Eqs.  (2)-(5))  are  transformed  into  the  frequency 
domain.  The  detailed  procedure  of  transformation  has  been 
explained  in  a  previous  study  [12].  In  contrast  to  modelling 
the  EIS  with  equivalent  circuit  models,  the  model  presented 
here  has  direct  physico-chemical  relevance,  since  it  does  not 
contain  empirical  elements.  Furthermore,  the  equations  and 
parameters  are  identical  to  that  for  steady-state  and  dynamic 
simulations. 


3.2.  Applied  method  of  parameter  identification 


Eq.  (1)  is  the  methanol  mass  balance  in  the  anode  compartment. 
Eqs.  (2)-(4)  are  the  mass  balances  of  methanol,  the  interme¬ 
diate  CO  and  the  intermediate  OH  in  the  anode  catalyst  layer. 
Eqs.  (5)  and  (6)  are  the  charge  balances  on  the  anode  and  the 
cathode  sides,  respectively.  Eqs.  (7)-(10)  display  the  reaction 
rates  at  the  anode  and  cathode.  Finally,  Eq.  (11)  describes  the 


To  determine  a  realistic  parameter  set  for  a  given  MEA,  the 
effects  of  the  various  reaction  and  transport  processes  inside  the 
MEA  are  separated  and  the  characteristic  parameters  are  iden¬ 
tified.  The  motivation  and  procedure  applied  is  as  follows.  A 
characteristic  feature  of  each  reaction,  charge  and  transport  pro¬ 
cess  is  its  specific  response  time  to  changes  in  the  state  of  the 
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Input:  Iqqh  Model  Output:  Ucen,  Z 


Fig.  2.  Principle  of  modelling  steady-state,  dynamic  and  impedance  measurements  using  one  equation  and  parameter  set. 


MEA.  Each  process  can  therefore  be  attributed  to  an  investiga¬ 
tion  method  which  covers  the  distinct  time  range  of  the  process. 
The  schematic  in  Fig.  3  illustrates  that  charge  processes  and 
fast  kinetics  influence  the  cell  performance  within  seconds  or 
milliseconds,  while  reactant/product  transport  or  temperature 
change  effects  are  observed  over  several  seconds  to  minutes. 
In  particular,  the  fast  electrochemical  processes  at  the  elec¬ 
trode  can  be  characterized  with  EIS.  In  previous  studies  of  the 
response  of  DMFCs  to  step  changes  in  current,  it  was  found 
that  a  combination  of  fast  kinetics  and  slow  kinetics  causes 
the  experimentally  observed  strong  voltage  response.  Transport 
processes  such  as  methanol  crossover  contribute  only  a  small 
additional  effect  after  the  main  dynamic  response  has  faded 
[9].  Finally,  steady-state  experiments  like  IV-measurements  and 
mass  balances  are  influenced  by  all  processes  in  their  steady- 
state.  These  experiments  are  especially  useful  to  determine  the 
transport  parameters,  since  transport  effects  occur  only  after  sev¬ 
eral  seconds.  They  are  barely  covered  by  the  other  investigation 
methods. 

For  this  research,  the  anode  kinetic  constants,  capaci¬ 
tances  and  surface  concentrations  were  approximated  by  fitting 
to  anode  impedance  spectra  and  current  step  responses.  Rei 
was  determined  from  the  anode  EIS  at  high  frequencies. 
GDF  mass  transport  constants  and  membrane  mass  trans¬ 
port  constants  were  fitted  to  polarization  and  mass  balance 
experiments,  respectively.  Finally,  the  cathode  reaction  kinetic 
constants  were  approximated  by  fitting  to  the  polarization 
curve. 


charge 

- 1 — ► 

fast  kinetics 

slow  kinetics 

species  transport 

temperature 

EIS 

steady  state 

current  step  response 

Fig.  3.  Velocity  of  processes  inside  the  DMFC  and  specific  time  range  covered 
by  experimental  investigation  methods. 


4.  Identification  of  anode  kinetic  parameters 

In  order  to  focus  on  the  influence  of  the  variation  of  the  dif¬ 
fusion  layer  and  membrane  on  a  MEA’s  performance,  one  MEA 
was  selected  and  the  parameters  related  to  its  catalyst  layer  were 
identified.  This  MEA  consists  of  Nation  115,  31BC  as  anode 
GDL,  the  modified  10DA  as  cathode  GDL  (see  Section  2.1),  the 
standard  cathode  catalyst  layer  (see  Section  2.1)  and  an  anode 
catalyst  layer  of  6  mg  cm-2  Pt-Ru  loading  with  a  Nafion/catalyst 
weight  ratio  of  0.07.  The  anode  catalyst  layer  differs  slightly 
from  that  used  when  conducting  the  diffusion  layer  and  mem¬ 
brane  transport  parameter  studies  (Section  5),  but  the  steady  state 
and  dynamic  behaviour  were  very  similar. 

Since  this  work  is  conducted  before  an  in-depth  study  of  the 
effect  of  membrane  and  diffusion  layers  on  MEA  behaviour,  the 
additional  set  of  equations  for  a  separate  description  for  MPL  and 
GDL  was  not  taken  into  account,  and  literature-based  transport 
parameters  were  used.  The  diffusion  coefficient  of  methanol  in 
Nation  115  (=6.26  x  10_10m2  s-1)  was  calculated  using  mea¬ 
surements  presented  by  Kallio  et  al.  [3].  The  effective  mass 
transport  coefficient  in  the  anode  GDL  (=1.154  x  10  5ms  l) 
was  approximated  by  combining  geometry  and  literature  values 
as  discussed  in  [9].  The  anode  flow  rate  was  10  x  10-6  m3  min-1 
and  the  active  area  was  10  x  10-4  m2. 

The  kinetic  parameters  were  fitted  iteratively  to  the 
impedance  spectra,  current  step  experiments  and  anode  and 
cell  polarization  experiments  as  motivated  in  Section  3.  The 
impedance  spectra  were  fitted  roughly  to  identify  a  suitable 
range  for  the  reaction  rate  constants  and  other  parameters.  It 
is  assumed  that  the  lumped  description  of  the  electrode  and  its 
structure  may  have  prevented  quantitative  reproduction  of  the 
experiments.  Nevertheless,  the  order  of  magnitude  and  ratio  of 
the  kinetic  and  charge  parameters  were  identified  (for  parameter 
sensitivity  see  [12]).  Simulation  of  the  current  step  responses 
showed  high  sensitivity  to  the  parameters  of  the  slow  reaction 
kinetics  (rAio,  ^A30,  geo,  cpt  and  cru),  as  predicted  in  Section 
3,  and  this  led  to  further  refinement  of  the  parameters. 

The  results  of  the  parameter  estimation  are  shown  in  compari¬ 
son  with  the  experimental  results  in  Figs.  4-6.  The  IV-curve  can 
be  reproduced  quantitatively  (Fig.  4).  Furthermore,  the  anode 
polarization  curve  was  very  close  to  the  experimental  one.  The 
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Fig.  4.  Steady-state  results:  experimental  (x)  and  simulated  (-)  steady-state  IV- 
curve;  experimental  anode  polarization  (O)  and  simulated  overpotentials  ( — ); 
surface  coverages  ( — ). 


experimental  data  points  of  the  latter  were  taken  during  the 
EIS  measurements.  The  simulation  result  predicts  the  following 
contributions  from  the  anode  and  the  cathode  sides: 

•  /cell  <  /cell, max :  similar  anode  and  cathode  loss 

•  /cell  — >  /cell, max'  anode-induced  limitation  (see  $oh  increase). 

Fig.  5  shows  that  the  model  also  quantitatively  predicts  the 
experimental  dynamic  behaviour  of  the  DMFC  to  various  step 
changes  in  current.  Since  the  step  response  is  governed  by 
anode  kinetics  [13]  and  even  the  steady-state  anode  potential 
is  reproduced  (Fig.  4),  a  reasonable  anode  parameterization  can 
be  assumed.  However,  as  explained  before,  the  same  parameter 
set  did  not  yield  a  quantitative  reproduction  of  the  experimental 
anode  EIS.  This  is  shown  in  Fig.  6.  While  the  shape  and  cur¬ 
rent  level  influence  are  qualitatively  reproduced,  the  simulated 
curves  are  shifted  to  lower  frequencies  (Bode  plot)  and  show 


Fig.  5.  Experimental  (x)  and  simulated  (-)  response  to  cell  current 
steps:  1800^  1000 Am-2  (top),  1800^  1400 Am-2  (middle)  and  3800^ 
3000  A  m-2  (bottom). 
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Fig.  6.  Experimental  (x)  and  simulated  (-)  electrochemical  anode  impedance 
spectra:  1000  Am-2  (black),  2000  Am-2  (grey)  and  3000  Am-2  (light  grey), 
(a)  Nyquist  plot  and  (b)  Bode  plot.  Impedance  spectra  were  corrected  for  ohmic 
drop. 


a  smaller  dependence  on  the  quasi- steady  state  cell  current  in 
comparison  with  the  experiment  (Nyquist  and  Bode  plot). 

5.  Membrane  and  anode  GDL  key  parameters 

5.7.  Membrane 

Determination  of  the  membrane  and  key  GDF  parameter  val¬ 
ues  was  done  simultaneously.  The  GDF  parameters  used  in  this 
section  for  membrane  characterization  are  those  presented  in 
Section  5.2  below. 

The  influence  of  the  membrane  on  DMFC  performance 
is  twofold;  ohmic  resistance  of  the  membrane  RQ\  causes 
ohmic  losses  in  the  DMFC  polarization  curve  and  furthermore, 
methanol  permeability  of  the  membrane  affects  the  polarization 
curve  and  Faradaic  efficiency.  The  methanol  permeability  of  a 
membrane  depends  on  the  membrane  thickness  dM  and  on  the 
diffusion  coefficient  of  methanol  in  the  membrane,  D^H3qH.  The 
contribution  of  electro-osmotic  drag  to  the  methanol  crossover 
flux  is  small  in  the  targeted  low  current  density  operating  range. 
Also,  the  electro-osmotic  drag  would  cause  a  positive  devia¬ 
tion  from  the  diffusion-induced  linear  decrease  of  crossover 
with  current  density.  Such  a  positive  deviation  has  not  been 
observed  in  the  literature  [10,11].  Contributions  by  electro- 
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osmotic  drag  are  therefore  not  accounted  for  by  the  model. 
Water  permeability  is  an  important  issue  [1,14]  but  it  is  not 
considered  in  this  state  of  the  ME  A  design  model.  The  three 
parameters  R& j,  dM  and  ^3qh  cover  the  two  main  influences 
on  the  membrane,  namely,  electrical  resistance  and  methanol 
crossover.  They  do  not  give  a  mechanistic  explanation  of  the 
observed  membrane  behaviour,  because  such  an  approach  is  suit¬ 
able  for  in-depth  material  specific  studies  but  not  for  setting  up 
a  database  consisting  of  arbitrary  membrane  materials,  as  tar¬ 
geted  in  this  work.  Instead,  the  user  of  the  DMFC  MEA  design 
model  should  be  familiar  with  all  the  parameters  and  be  able  to 
measure  or  estimate  them  easily.  The  membranes  investigated 
and  described  here  are  a  commercial  Nation  115  membrane  and 
three  hydrocarbon-based  membranes;  they  are  listed  in  Table  1. 


Table  1 

Key  membrane  parameters  of  the  DMFC  MEA  model 


Membrane 

dM  (wet) 

(x 10-6  m) 

Re l  (x  10-4  £2  m2) 

Dm 

^ch3oh 

(xlO-10m2  s-1) 

N115 

142 

0.133 

5 

Ml 

43.2 

0.065 

1.45 

M2 

39 

0.18 

0.8 

M3 

41 

0.13 

1 

For  all  membranes,  the  ohmic  resistance  of  MEAs  containing  the 
respective  membrane  was  measured  by  impedance  spectroscopy 
at  343  K.  It  should  be  noted  that  these  resistance  values  also 
contain  a  small  contribution  from  the  electronic  resistances  of 


MEA  [Ml,  AD:  31 BC,  CD:  10DA+MPL].  FA=3  ml/min,  1M  CH3OH 

—  sim.cell  voltage 

- sim.  anode  overpotential 

O  exp.  cell  voltage 
*•'  exp.  cell  voltage 


Current  density,/  [Am  ] 


MEA  [M2,  AD:  31  BC,  CD:  10DA+MPLJ.  FA  =3  ml/min,  1M  CH3OH 

—  sim.cell  voltage 

- sim.  anode  overpotential 

x  exp.  cell  voltage 
0  exp.  anode  overpotential 


Current  density, /'ce||  [A  m"  ] 


MEA  [M3,  AD:  35BC,  CD:  35DC] 


Current  density,/  [A  m'2]  Current  density,/  [A  m*2] 


M£A  [M3,  AD:  35BC,  CD:  35BC],  FA  =  1.7ml/min,  1M  CH3OH 


Current  density, /'ce||  [A  m’2] 


Fig.  7.  Steady-state  characterization  of  MEAs  with  membranes  Ml  (a  and  b),  M2  (c  and  d),  M3  (e  and  f).  Left:  polarization  curve;  anode  overpotential.  Right: 
equivalent  current  density  loss  by  methanol  crossover.  MEA  composition  (AD:  anode  diffusion  layer,  CD:  cathode  diffusion  layer)  and  operating  conditions  are  given 
inside  plots. 
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the  respective  MEA.  The  membrane  thickness  was  measured 
experimentally  after  exposure  to  1  M  methanol  solution.  The  dif¬ 
fusion  coefficient  was  adjusted  so  that  the  simulated  methanol 
crossover  flux  of  the  respective  MEA  matches  the  experimental 
one.  A  comparison  of  the  experimental  and  simulated  results  is 
shown  in  Fig.  7  for  ME  As  containing  the  membranes  Ml,  M2 
and  M3  and  in  Fig.  8  for  various  ME  As  containing  Nation  115. 
The  exact  values  of  the  key  parameters  for  each  investigated 
membrane  are  listed  in  Table  1 . 

5.2.  Anode  diffusion  layer 

Determination  of  the  key  parameters  for  the  diffusion  layer 
is  more  complex  than  that  for  the  membrane.  The  effective 
methanol  mass  transport  coefficient  varies  with  PTFE  content, 


porosity  (or  alternatively  areal  weight),  thickness,  and  type  of 
micro-porous  layer  (MPL)  [2,15,16].  The  existence  or  absence 
of  an  MPL  significantly  influences  the  mass  transport  to  the 
anode  catalyst  layer,  as  shown  in  Fig.  8b  and  d;  the  MEA  without 
MPL  shows  significantly  higher  methanol  crossover  (Fig.  8b) 
than  the  MEA  with  MPL  (Fig.  8d).  Influential  MPL  parame¬ 
ters  are  thickness,  loading  level  and  PTFE  content.  For  a  given 
MEA,  the  values  for  all  six  parameters  are  known  or  can  alter¬ 
natively  be  obtained  by  experiment.  The  PTFE  content  of  GDL 
and  MPL,  GDL  porosity,  GDL/MPL  combined  thickness  and 
MPL  loading  level  are  provided  by  the  manufacturer.  As  an 
alternative  to  the  GDL  porosity,  the  GDL’s  areal  weight  can  be 
measured.  MPL  and  GDL  thicknesses  are  obtained  using  a  light 
microscope  as  shown  in  Fig.  9  (video  microscope  system  Some- 
tech  SV  35).  Microscopic  pictures  of  the  GDL/MPL  3 1DC  show 


(a)  MEA  [Nafionl  1 5,  AD:  31  BA,  CD:  31 BC].  F.  =3  ml/min 

0.7j- 

- sim.  cell  voltage.  1M  CH3OH 

- sim.  anode  overpotential,  1 M  C 

- sim.  cell  voltage.  0.5M  CH3OH 

•x  exp.  cell  voltage.  1M  CH3OH 

exp.  anode  overpotential.  1M  CH^OH 
4-  exp.  cell  voltage.  0.5M  CH3OH 


6000  8000 

Current  density,  /'ce||  [A  m'2] 


(b)  MEA  [Nafionl  1 5,  AD:  31  BA,  CD:  31  BC],  FA  =3  ml/min,  1 M  CH3OH 


(C)  MEA  [Nafionl  1 5,  AD:  31  BC,  CD:  31  BC], 


cell  voltage,  1M  CH3OH 


exp.  cell  voltage,  1M  CH3OH 


+  exp.  cell  voltage.  0.75M  CH3OH 
O  exp.  cell  voltage.  0.5M  CH3OH 

- sim.  cell  voltage.  0.5WI  CH3OH 

- sim.  cell  voltage.  0.75M  CH3OH 


6000 


(e) 


0.7 


Current  density,  /'ce||  [A  m  ] 

MEA  [Nil 5,  AD:  31  DC,  CD:  31  BC],  FA=3  ml/min 
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—  sim.  cell  voltage,  1 M  CH3OH 

- sim.  anode  overpotential,  1M  CH^OH 

x  exp.  cell  voltage,  1M  CH3OH 
O  exp.  anode  overpotential,  1 M  CH  OH 


1000 


2000 


3000 


4000 


Current  density,  /'ce||  [A  m  J 


(f)  MEA  [N115,  AD:  31  DC,  CD:  31  BC],  FA  =3  ml/min,  1 M  CH3OH 


Fig.  8.  Steady-state  characterization  of  MEAs  with  Nation  115  and  various  anode  diffusion  layers:  3 1BA  (a  and  b),  3 1BC  (c  and  d),  3 1DC  (e  and  f).  Left:  polarization 
curve;  anode  overpotential.  Right:  equivalent  current  density  of  methanol  crossover.  MEA  composition  (AD:  anode  diffusion  layer,  CD:  cathode  diffusion  layer)  and 
operating  conditions  are  given  inside  plots. 
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Table  2 


Key  parameters  for  diffusion  layer  and  micro-porous  layer  (MPL) 


Anode  GDL  and  MPL 

dAD  (excl.  MPL) 

(x 10-6  m) 

^"PTFE  (®XCl. 

MPL)  (%) 

AWad  (excl. 

MPL)  (kgm-2) 

eAD  (excl.  MPL) 

jAMPI. 

(x 10-6  m) 

*ptfeL 

AMpL 

carbon 

(mg  cm-2) 

31BA 

310 

5 

0.065 

- 

- 

- 

- 

31BC 

270 

5 

0.065 

- 

70 

23 

4 

31DC 

270 

20 

0.0743 

- 

70 

23 

4 

35BA 

295 

5 

0.0548 

- 

- 

- 

- 

35BC 

270 

5 

0.0548 

- 

70 

23 

4 

35DC 

239 

20 

0.0626 

- 

70 

23 

4 

25BC 

185 

5 

- 

0.88 

70 

23 

4 

10DA  +  SDI  MPL: 

1.3  mg  cm-2 

380 

20 

0.085 

— 

40 

50 

1.3 

the  GDL  and  MPL  thickness  (Fig.  9c)  and  highlight  the  strong 
structural  differences  between  the  porous  GDL  (Fig.  9a)  and  the 
dense  MPL  (Fig.  9b).  An  overview  of  the  exact  values  of  the 
key  parameters  for  several  GDL/MPLs  is  given  in  Table  2.  In 
order  to  keep  the  DMFC  MEA  design  model  usable  by  many 
users  as  well  as  extendable  for  new  materials,  modelling  of  the 
diffusion  layer  is  undertaken  with  only  the  above-mentioned  key 
parameters  as  input  parameters. 

The  diffusion  layer  mass  transfer  coefficient  is  calculated  as 
follows: 

gAD 

&AD  =  DCH30H-^  (15) 

4P 

If  the  diffusion  layer  porosity  £AD  is  not  known,  it  is  calculated 
by  the  experimentally  measurable  areal  weight  AWad,  the  PTFE 
content  of  the  diffusion  layer  xFFFE  and  the  bulk  densities  of  the 
diffusion  layer  components  carbon,  pcarbon  =  1.7  x  103  kgm-3 
and  PTFE  pptfe  =  2.19  x  103  kgm“3  [17]: 


AWAD  (  rAD  1  —  rAD  \ 

AVV  /  XPTFE  _j_  1  APTFE  j 

dAD  y  PPTFE  Pcarbon  J 


(16) 


Pcarbon  was  determined  by  applying  Eq.  (16)  to  a  set  of  GDLs 
with  0  wt.%  PTFE  content  (10AA,  3 1 AA  and  30 AA)  and  match¬ 
ing  the  given  values  for  AWad,  dAD  and  £AD. 

An  effective  diffusion  layer  thickness  dA®,  which  is  based  on 
the  original  key  parameter  dAD ,  is  introduced;  dA^  accounts  for 
MEA  compression  and  for  the  pressure-driven  convective  flow 
between  adjacent  anode  channels  through  the  GDL.  Accord¬ 
ing  to  Ref.  [6],  the  pressure-driven  convective  flow  through  the 
GDL  close  to  the  anode  flow-field  ribs  significantly  increases 
methanol  mass  transport  inside  the  diffusion  layer.  Also,  calcula¬ 
tions  with  the  model  presented  above  indicate  that  an  additional 
transport  phenomenon  besides  diffusion  has  to  be  present.  When 
maximizing  the  DMFC  model’s  GDL  diffusion  transport  by 
setting  the  GDL  porosity  to  1  and  keeping  the  original  GDL 
thickness,  the  model  predicts  that  the  anode  methanol  concentra¬ 
tion  reaches  0  mol  l- 1  at  current  densities  which  are  significantly 
lower  than  the  experimentally  obtained  limiting  current  densi¬ 
ties.  To  account  for  the  pressure-induced  convective  transport,  a 
hydrodynamic  layer  inside  the  GDL,  which  is  governed  by  con¬ 
vective  flow,  is  assumed.  The  layer  decreases  the  real  diffusion 


layer  thickness  to  an  effective  value  for  diffusive  mass  trans¬ 
port,  d^f.  For  a  50%  compressed  diffusion  layer  with  5  wt.% 
PTFE  and  thickness  dAD  >280  x  10-6  m,  a  constant  thickness 
of  the  hydrodynamic  layer  (dA^  =  60  x  10-6  m)  is  identified, 
as  explained  in  the  discussion  of  Eq.  (17).  For  thicknesses 
dAD<  280  x  10  6  m,  is  assumed  to  decrease  linearly  to 
zero  until  dAD  =  0  m.  In  addition,  the  hydrodynamic  layer  inside 
the  GDL  depends  on  the  PTFE  content  of  the  GDL.  A  high 
PTFE  content  increases  GDL  hydrophobicity  and  mechanically 
hinders  convective  flow  due  to  formation  of  thin  PTFE  planes 
inside  the  GDL.  Both  effects  should  decrease  the  mass  transport 
through  the  hydrodynamic  layer.  A  linear  influence  of  PTFE 
content  on  the  hydrodynamic  layer  thickness  is  assumed.  The 
effective  diffusion  layer  thickness  is  then  calculated  as  follows: 


./AD 

\ff 


0.5  d 


AD 


./AD 

3hyd 


aptfe  0-05  \ 
0.2-0.05  J 


(17) 


Here  0.2  and  0.05  are  the  PTFE  contents  (wt.%)  in  the  anode 
GDLs  of  the  MEAs  presented  in  Fig.  8.  Validation  of  Eq.  (17) 
was  done  using  the  experimental  polarization  curves  presented  in 
Fig.  8c  and  e.  The  experimental  limiting  current  densities  shown 
in  Fig.  8c  and  e  are  strongly  dependent  on  the  PTFE  content  and 
are  reproduced  by  the  model  for  several  concentrations.  The 
dependence  of  methanol  loss  on  concentration  is  also  correctly 
predicted,  but  the  experiment  shows  a  lower  actual  effect  of 
PTFE  content  on  methanol  crossover. 

The  MPL  mass  transfer  coefficient  is  calculated  as  follows: 
gAMPL 

k  =  Dch3oh  tampljampl 
with  the  assumption 


cAMPL  .  _ 

b_ _  _  /  AMPLx1-5 

t-AMPL  } 


(19) 


Here  £AMPL  is  the  MPL  porosity,  rAMPL  the  MPL  tortuosity, 
and  dAMPL  is  the  MPL  thickness.  Eq.  (19)  is  used  for  micro- 
porous  layers  such  as  the  catalyst  layer  or  MPL.  MPL  porosity 
is  calculated  using  the  MPL  key  parameters: 


gAMPL 


jAMPL  _  (w;AtMfPeL/pptfe)  _  (u;AXP  eff/Pcarbon) 
jAMPL 


(20) 
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Fig.  9.  Micrographs  of  GDL/MPL  31  DC.  (a)  View  of  uncoated  side  of 
GDL/MPL,  (b)  view  of  MPL-coated  side  of  GDL/MPL  and  (c)  cross-sectional 
view  of  GDL/MPL  with  MPL  on  upper  side. 


where  w^oneff  the  effective  carbon  loading  of  the  MPL 
and  Wp|p£L  the  PTFE  loading  calculated  from  the  PTFE  mass 
fraction  : 


w 


AMPL 

PTFE 


AMPL 
^carbon,  eff 


AMPL 
VpTFE 
i  AMPL 
1  XPTFE 


(21) 


Modelling  of  the  MPL  effect  by  using  the  original  MPL 
key  parameters  requires  some  adjustments.  An  equivalent  pro¬ 
cedure  had  been  introduced  when  modelling  the  GDL  effect 
on  DMFC  performance  with  GDL  key  parameters.  Using  the 


original  MPL  loading  and  the  MPL  thickness  determined  by 
microscopy,  the  model  predicts  insufficient  methanol  diffusion. 
As  a  result,  the  simulated  limiting  current  densities  were  sig¬ 
nificantly  lower  than  the  experimental  ones.  Discussion  with 
MPL  producers  and  analysis  of  micrographs  indicate  that  not  all 
the  applied  MPL  material  is  located  in  the  MPL  layer.  During 
the  MPL  manufacturing  process,  the  macro-pore  sized  diffusion 
layer  permits  the  transport  of  a  certain  amount  of  MPL  mate¬ 
rial  deep  into  the  diffusion  layer;  some  was  even  reported  to 
have  completely  passed  the  diffusion  layer.  Comparing  experi¬ 
mental  and  simulated  results  of  MEAs  with  and  without  MPL 
(Fig.  8a-d)  indicates  that  the  effective  carbon  loading  may  be  60 
to  100%  of  the  original  carbon  loading.  This  value  is  assumed 
to  depend  on  PTFE  loading  in  the  GDL  since  this  will  affect 
penetration  of  the  MPL  material  into  the  GDL.  The  following 
equation  is  proposed  for  the  effective  carbon  loading: 

„  AMPL  _  „,,AMPL  (a/;  ,  a/I  XPTFE  —  0'^  \ 

^carbon,  eff  ^carbon  10.6  +  0.4  Q2_QQ5  j  @2) 

Fig.  8a  and  b  shows  the  performance  of  a  ME  A  without 
MPL  and  Fig.  8c  and  d  shows  that  of  a  MEA  with  an  MPL. 
The  predicted  limiting  current  densities  are  close  to  those  of 
the  experiments.  The  concentration  dependence  is  also  correctly 
reproduced.  Probably  due  to  gaseous  carbon  dioxide  blocking 
the  methanol  transport  at  high  current  densities,  however,  the 
experimental  cell  voltage  is  lower  than  that  predicted;  this  holds 
also  for  the  MEA  with  20%  PTFE  content  shown  in  Fig.  8e 
and  f.  Finally,  comparing  the  simulation  results  of  the  diffusion 
layer  and  the  influence  of  the  membrane  on  the  corresponding 
experimental  results,  as  shown  in  Figs.  7  and  8,  it  can  be  con¬ 
cluded  that  the  model,  using  a  basic,  experimentally  oriented  set 
of  key  parameters  for  the  diffusion  layer  and  membrane,  is  able 
to  predict  the  performance  of  a  variety  of  MEAs. 

6.  Case  studies  for  membrane  and  diffusion  layer  design 

Since  high  Faradaic  efficiency  is  a  crucial  issue  in  design¬ 
ing  DMFC  systems,  MEAs  with  minimized  methanol  crossover 
are  highly  advantageous.  To  quantify  the  effect  of  membranes 
on  DMFC  performance,  a  case  study  using  existing  mem¬ 
branes  is  presented  in  the  following.  Fig.  10  gives  a  comparison 
between  the  performance  of  an  MEA  containing  one  M3  mem¬ 
brane  and  an  MEA  containing  a  stack  of  two  M3  membranes. 
As  can  be  seen  from  Table  1,  M3  already  has  a  significantly 
lower  methanol  diffusion  coefficient  (D^H3qH)  than  Nation  115. 
This  positive  material  property  is  neutralized,  however,  by  a 
much  lower  membrane  thickness  compared  with  Nation  115, 
which  leads  to  similar  mass  transfer  coefficients  (^ch3oh/^M) 
for  both  membranes.  Using  two  membranes,  or  equivalently 
a  membrane  M3  with  a  double  thickness,  should  decrease 
methanol  crossover  but  will  also  increase  ohmic  resistance 
twofold.  The  simulations  were  done  at  a  stoichiometric  methanol 
feed  of  2.5  and  a  methanol  inlet  concentration  of  1  M.  Fig.  10a 
shows  that  the  MEA  with  two  membranes  has  a  slightly  lower 
cell  current  density  for  DMFC  operation  at  0.45  V.  The  dif¬ 
ference  in  electrochemical  performance  between  both  MEAs 
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Fig.  10.  Performance  of  MEA  with  one  M3  membrane  (dashed  line)  and  with 
two  M3  membranes  (solid  line)  at  anode  stoichiometry  2.5.  (a)  Cell  voltage 
(black)  and  methanol  loss  (grey),  (b)  Methanol  concentration  profile  inside 
MEA.  MEA  composition  (AD:  anode  diffusion  layer)  and  operating  conditions 
are  given  inside  plots. 

increases  to  significant  amounts  for  high  current  densities. 
While  the  MEA  with  two  membranes  shows  a  decrease  in 
electrochemical  performance  as  discussed,  a  30%  reduction  in 
methanol  crossover  flux,  is  achieved.  The  proportional  decrease 
in  methanol  crossover  flux,  which  is  achieved  when  using  the 
MEA  with  two  membranes,  is  similar  for  all  current  densities. 
By  contrast,  the  difference  between  the  absolute  amounts  of 
methanol  lost  by  crossover  becomes  insignificant  at  high  cur¬ 
rent  densities.  Fig.  10b  shows  the  influence  of  the  different 
membranes  on  the  methanol  concentration  profile  in  the  MEA. 
The  MEA  with  two  membranes  significantly  increases  the  mean 
methanol  concentration  level  inside  the  anode  catalyst  layer 
and  leads  to  a  slight  increase  in  methanol  concentration  even 
inside  the  anode  channel.  The  low  electrochemical  performance 
at  high  current  densities  and  the  rather  slight  improvement  in  the 
corresponding  crossover  make  the  MEA  with  two  membranes 
advantageous  only  for  operation  at  low  current  densities  and 
high  cell  voltages. 

High  anode  flow  rates  correspond  to  high  water  loss  and 
high  electricity  consumption  by  the  feed  pump;  low  anode  flow 
rates  may  therefore  be  favourable  for  DMFC  systems.  If  lower 
anode  flow  rates  are  used  and  the  same  anode  feed  stoichiom¬ 


MEA  [2XM3,  AD:  35BC],  JL  =  2.5 


Fig.  1 1 .  Performance  of  MEA  with  two  M3  membranes  for  1  M  CH3OH  (solid 
line)  and  for  3  M  CH3OH  (dashed  line)  inlet  concentration  at  anode  stoichiom¬ 
etry  2.5.  (a)  Cell  voltage  and  methanol  loss,  (b)  Methanol  concentration  profile 
inside  MEA. 


etry  is  maintained,  the  inlet  methanol  concentration  has  to  be 
increased.  Fig.  1 1  shows  a  comparison  of  MEA  performances 
at  methanol  inlet  concentrations  of  1  and  3  M;  the  MEA  is  the 
one  discussed  previously,  with  two  M3  membranes  and  which 
shows  improved  performance  in  the  high  voltage  region.  Using 
inlet  methanol  concentrations  of  3  M  instead  of  1  M  leads  to 
slightly  increased  electrochemical  performances,  but  to  greatly 
increased  methanol  crossover  (Fig.  11a);  methanol  loss  increases 
by  more  than  twofold  when  using  a  3  M  methanol  solution  and 
is  even  higher  than  that  for  the  MEA  with  a  single  M3  at  1  M 
(Fig.  10).  Fig.  lib  illustrates  that  the  level  of  methanol  concen¬ 
tration  inside  the  MEA  for  3  M  inlet  concentrations  is  between 
1.5  and  1  M,  in  other  words  more  than  half  of  the  original  inlet 
concentration.  The  main  reason  for  the  increased  methanol  loss 
for  the  same  stoichiometric  flow  rate  is  the  transport  of  methanol 
inside  the  MEA  driven  by  the  concentration  difference.  It  can 
be  concluded  that  increasing  membrane  thickness  is  beneficial 
for  DMFC  operation  at  low  current  densities,  while  an  increase 
of  the  inlet  methanol  concentration  has  a  strongly  negative 
impact  on  the  Faradaic  efficiency  when  operating  at  low  current 
densities.  This  holds  even  for  membranes  with  low  methanol 
permeability. 
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Way  x  in  MEA  [pm*1]  Way  x  in  MEA  [pm*1] 

Fig.  12.  Performance  of  MEA  with  modified  anode  MPL  and  two  M3 
membranes  for  1M  methanol  and  anode  stoichiometry  2.5:  35BC  (solid 
line);  35BC  +  additional  MPL  (2  mg  cm-2  carbon  loading  +  35  [xm  thickness 
increase),  (a)  Cell  voltage  and  methanol  loss,  (b)  Methanol  concentration  profile 
inside  MEA. 

A  final  case  study  should  show  the  influence  of  increased 
MPL  loading  on  performance.  In  order  to  increase  further 
the  performance  of  the  MEA  with  two  membranes  (see 
Figs.  10  and  11),  the  MPL  loading  of  the  anode  diffusion  layer 
is  increased  1.5  times.  It  is  assumed  that  the  MPL  thickness 
also  increases  1.5  times.  The  respective  results  are  presented  in 
Fig.  12.  As  shown  in  Fig.  12a,  the  MEA  with  the  additional  MPL 
loading  displays  the  expected  drop  in  limiting  current  density, 
whereas  the  electrochemical  performance  is  essentially  identi¬ 
cal  at  0.45  V  and  above.  The  methanol  loss  decreases  for  all 
current  densities,  but  notably  more  at  higher  current  densities. 
The  methanol  concentration  profile  at  0.45  V  (Fig.  12b)  illus¬ 
trates  the  effect  of  the  additional  MPL  loading.  Due  to  a  steeper 
concentration  drop  inside  the  MPL,  the  mean  methanol  concen¬ 
tration  inside  the  anode  catalyst  layer  decreases  and  causes  less 
methanol  crossover.  It  can  be  concluded  that,  when  operating  at 
0.45  V,  the  MEA  with  additional  MPL  loading  may  be  beneficial 
since  it  shows  improved  methanol  efficiency  while  maintaining 
the  same  electrochemical  performance. 

7.  Conclusions 

This  research  proposes  a  DMFC  MEA  design  model.  The 
comprehensive  model  presented  here  is  designed  to  be  easily 
extendable  so  that  additional  phenomena  can  be  incorporated 


at  a  later  stage.  Special  attention  is  given  to  the  selection  of 
model  input  parameters  which  are  easily  obtained  experimen¬ 
tally;  this  aspect  facilitates  building-up  and  future  enlargement 
of  the  material  database  as  well  as  extending  the  usability  of  the 
model. 

For  parameter  identification,  MEAs  are  investigated  using 
different  experimental  methods  such  as  polarization,  current- 
steps  and  EIS.  Since  the  experiments  cover  different  response 
times  of  the  DMFC,  processes  in  the  DMFC  such  as  the 
transport,  reaction  and  charge  processes  can  be  separately  inves¬ 
tigated.  The  main  reaction  kinetic  parameters  are  identified  at 
low  current  densities  using  anode  EIS  and  current- step  exper¬ 
iments.  By  systematic  variation  of  the  properties  of  the  single 
layers  inside  the  MEA,  it  is  possible  to  identify  key  parame¬ 
ters  for  reproducing  the  influence  of  each  respective  layer  on 
MEA  performance.  In  this  research,  the  basic  structure  of  an 
extendable  MEA  design  model  is  built  up  and  variations  of 
the  properties  of  the  anode  diffusion  layer  and  membrane  are 
integrated.  The  model  is  presented  with  a  basic,  experimen¬ 
tally  oriented  set  of  key  parameters  and  is  able  to  predict  the 
performance  of  a  variety  of  MEAs. 

Case  studies  for  membrane  and  diffusion  layer  designs 
demonstrate  the  potential  usefulness  of  the  model  for  MEA 
design.  The  model  predicts  that  increasing  the  membrane  thick¬ 
ness  decreases  methanol  crossover,  especially  at  low  current 
densities,  while  additional  MPL  loading  decreases  methanol 
crossover,  especially  at  high  current  densities.  For  DMFC  oper¬ 
ation  around  0.45  V,  however,  both  modifications  are  beneficial 
for  increasing  the  Faradaic  efficiency  while  changes  in  the  elec¬ 
trochemical  performance  are  negligible.  In  contrast  to  these 
findings,  operation  at  higher  methanol  concentrations  at  cell 
voltages  of  0.45  V  results  in  significantly  increased  methanol 
crossover.  This  increase  is  predicted  to  even  neutralize  the  ben¬ 
efits  obtained  by  improvements  of  the  membrane  or  diffusion 
layer.  For  operation  at  high  voltages,  it  is  therefore  recom¬ 
mended  to  adjust  the  inlet  methanol  concentration  to  around 
1  M  or  less. 
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